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We consider the three-dimensional electron gas confined by a strictly two-dimensional homoge- 
neous positive charge density at 2 = 0. Within the Hartree-Fock approximation, we study the mode 
structure in the confined direction in the metallic regime. We find, that for r s < 1.3 (r s < 2.5) the 
unpolarized (polarized) electron gas starts to populate also the first excited state in the z-direction. 
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I. INTRODUCTION 

The two-dimensional homogeneous electron gas (2DEG) is one of the most simple and thus widely used model 
to study electronic correlations in two dimensions 1 ' 3 . Experimentally, two-dimensional electronic systems have been 
realized using heterostructures, e.g. semiconductor-insulator interfaces, where layers of electrons are tightly confined 
in one spatial dimension (z) by strong surface electric fields, and the discreteness of the quantized energy levels in 
z becomes important!^. However, since electronic wave functions and electromagnetic fields spatially extend in the 
z-direction, theoretical predictions for the 2DEG must be modified before a quantitative comparison is possible^. 

In order to study general effects due to the interplay of electron-electron interactions and correlations with the finite 
extension of the electronic density in the z-direction, we introduce the model of a quasi-two-dimensional electron gas 
(Q2DEG). This model provides a simple and natural extension of the 2DEG which contains essential features of more 
sophisticated microscopic descriptions of heterostructures^. Frequently, experiments are modeled with additional 
parameters to account for the finite thickness. In general, these parameters should not be considered as independent 
of the density due to charge neutrality. 

Similar to the electron gas in two and three dimensions, we consider a jellium of electrons in a positive charged 
background insuring total charge neutrality. Whereas the electrons are treated fully three-dimensional, the background 
charges remain strictly two-dimensional, described by a homogeneous charge density, a®, in the plane z = 0. For 
vanishing total charge of the system, the electrons are confined around the plane z — 0. Similar to the 2DEG, we 
introduce the dimensionless parameter r s = l/(asy / 7nTo) where as — h 2 / {m e e 2 ) is the Bohr radius, m e the mass 
and (— e) the charge of the electron. At zero temperature, the system is fully described by the value of r s which 
characterizes the effective two-dimensional density of the electrons. 

In this paper, we study the Q2DEG in the metallic density region (0.5 < r s < 5) in the Hartree-Fock approximation. 
In particular we determine the spatial density distribution of the electrons in the z-direction, and the possible transition 
between the occupation of a single confined mode to the occupation of two or more excited modes, or subbands, in z. 
We show that for r s — > the energy per particle can be written as: 

E m (c a , r.) = § £ cl + E ^fA + X ^ys) + Cm(caj rJ (1) 



where m is the number of occupied modes in the z-direction and £ m and X m are smooth functions of r s , determined 
within Hartree-Fock (HF), K p is a constant for fixed spin polarization, p, and c a are the concentrations of electrons in 
each mode. The correlation energy beyond Hartree-Fock, C m , is estimated within density-functional theory. At fixed 
density (fixed r s ), we determine the ground state for given concentrations, c a , and, finally, minimize with respect 
to the concentrations to obtain -E , m (c ajm i n , r s ). The main goal of this paper is to determine the density where the 
two-mode solution (section [v]) becomes energetically favorable compared to the single mode solution. 

The paper is organized as follows. Section [H] introduces the model Hamiltonian of the Q2DEG and discusses the 
technical problems related to the thermodynamic limit and the long range behavior of the Coulomb 1/r-potential in 
the potential energy. In Section |III| we use the HF approximation to simplify the many-body problem, and discuss 
the general structure of the ground state energy in the high density limit, r s — > 0. In the following sections, Section |TV| 
and SectiorjVj we discuss the single mode and two mode solution of the HF approximation. For both cases, we first 
start discussing the Hartree-approximation, where we have found analytical solutions for the resulting non-linear 
Schrodinger equation. These solutions serve to obtain a first estimate for the Hartree and exchange contribution 
to the energy, £^ and X^, respectively. We will show later, that the numerical minimization of the full HF-energy 
introduces only minor corrections. Finally, we briefly discuss correlation effects beyond HF within the local density 



approximation using density functional theory (Section VI I 



II. QUASI-TWO-DIMENSIONAL ELECTRON GAS MODEL 

Let us consider N electrons interacting with a homogeneous positive charged, strictly two-dimensional plane at 
z = and area S — L 2 . Assuming a charge-neutral system, the background surface density writes cro = N/S. The 
iV-body Hamiltonian is given by 

N h 2 

i— 1 

where Vn is the total potential energy of all charges. 

It is well known that the Coulomb potential poses difficulties in the definition of the potential energy in the 
thermodynamic limit due to the non integrability at infinity. The local singularity of the Coulomb potential near the 
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origin is a classical problem of self-adjointness and here we only focus on the definition of the potential with periodic 
boundary conditions. 

Let A denotes the two-dimensional lattice in ]R 3 generated by the vectors (I/,0,0) and (0,L, 0). For a regular 
integrable interaction v, we formally define the total periodic potential as: 



v N = v ee + v eb + v bb 
v K , 



i tGA,t^0 



dr dr v(r — r ) 



N 2 

.'sx 



(3) 
(4) 

(5) 

(6) 



where the index e holds for the electrons and b holds for the positive background. The last term in Eq|4] is the 
interaction of an electron with all its periodized images. As soon as the interaction v is regular and integrable, we 
can rewrite the potential energy as: 



V N =, 



with 



v(R) = y2 (v(R + t) - / rfr + t + r) 



«i(z) = / drv((r, z)) - v((r,0)) 



(7) 



and C„ is the Madelung energy of electrons on the lattice A in a homogeneous background 

[ V { T ) ~ ~q I drv(r + r)j- 



C v = ^2 (v(t)- Jjirv(r + T)J - jj I <Ir r( r) 



(8) 
(9) 

(10) 



Let us notice that the Fourier transform v(K) of v is directly related to the Fourier transform, v(K), of v. As can be 
directly verified, we have v(K) = v(K), except that v(K) = for k x = k y — 0. With this new definition (Eq{7]), we 
only need that v, V\ and C v are well defined. That is : 



V v(R + t)-\ [drv(R + 



< +oo 



dr \v((r, z)) - v((r, 0))| < +oo 



(11) 
(12) 



These conditions are fulfilled by the Coulomb potential Vq{R) = e 2 /R, except at R = as mentioned above. 
Furthermore, we have v\(z) = —2ire 2 \z\ and v(K) = Aire 2 /(k 2 + k 2 + k 2 ) for k x ^ or k y ^ and otherwise, 
and the periodic potential energy, Eq. finally writes 



v q2D {Ri 3 ) - 



27T \Zi - Zj\ 



dk z 47r 
2^ k 2 + k 2 ' 



S 



AKR 



27re 2 CT ^2 \ z i 



N 



C v 



(13) 
(14) 



where K = (k, k z ) and e lk ' T = 1 for r G A. 
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III. HARTREE-FOCK APPROXIMATION 



Within the HF approximation we minimize the ground state energy per particle, E, with respect to variations of 
the many-body wave-function, = det |{^if}| det ({^i^H, in the subspace of single Slater determinants 

In the following we assume that {ll^} (& = t>4-) are normalized, orthogonal single particle wave-functions, and we 
obtain for the total energy per particle: 



— — / dRdE! n e (R)v q2D {R- R')n e (R!) 
N 2 J SxR 

^yE / dRdR'^UR)^ 1 AR)v q2D {R-R')^,AR l )^ ] AR l ) 

i,j,a JSxR 

1 f ire 2 f 

— dRn e (R)2ire 2 a \z\- — - / dRdR' n e {R)\z - z'\n e {R!) 

^ JSxR iVi JSxR 



'Sxl 

ire 2 \ - 
"iVS ^ 



/ dRdR , yUR)v j *(R)\z-A*i*(R'W j A R ') + C T ( 16 ) 

where we have defined the total electronic density by 

n e (R) = J2\*i*(R)f (17) 

icr 

In this paper we are interested in a quasi-two-dimensional regime where we expect that the electrons populate a 
finite number, m, of discrete modes in the z-direction, whereas the density of states is continuous in the plane at 
constant z, in the thermodynamic limit. Each single-body wave function is then taken as a product of a plane 
wave 4>u in the plane z = and a wave function ip a in the z-direction where a labels the mode. Let N a be the number 
of electrons in the mode a and c a — N a /N, N — ^a- Accounting for the spin polarization p of the electrons, we 
have N a = iV a ^ + N a , and c a = c a ^ + c a ^. In the following, we restrict the discussion to i) the fully polarized gas 
(p = 1) where c a — c a ^ and c j, — 0, and ii) the unpolarized electron gas (p — 0) with c a ^ = c a ^ = c a /2 (unpolarized 
in each mode a). We further assume that the wave functions do not depend on spin: ip atT = %j) a . 

In each mode a, a, all transverse plane waves are occupied up to k F a — ^Jc aa k F with k F a B = 2/r s , and, in the 
thermodynamic limit, all summations over transverse states are replaced by integrals inside the Fermi surfaces 

5>E E -Ejr/ f*. as) 

aa \k\<k Fa<T aa %K F J \k\<kr a „ 



Further, for N — >• oo, the last line of Eq. (16) vanishes. 

It is instructive to regroup the different contributions to the total energy (in Hartree) as follows 

E[c a , M = + + *k»*» r '] (19) 

r s r s ' r s 

where the first term is the in-plane, strictly two-dimensional, kinetic energy with Kq = 1/2 for the unpolarized and 
K\ = 1 for the polarized electron gas. In order to separate the explicit r s -dependency in the following two terms, we 
introduce u = r^ kpZ together with the normalization J R du\tp ar:r (u)\ 2 — 1 of the confined modes. All contributions 
independent of the in-plane modes are contained in £ 

£[c a ,^a] = ~2^c Q f du^„(u)VC(u)+ / du p(u)^ + f du p{u)v p {u) (20) 

a JR JR L JR 

where the electrostatic potential, v p , is determined by the one-dimensional Poisson equation 

v'p(u) = 6(u) - P (u) (21) 
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from the total electronic density distribution p(u) = J2 acr c aa \ ipaa | 2 = J2 a c a \ip a \ 2 and the positive background charges 
at z = 0. Using v p {oo) — v'Joo) — 0, we have 

v Ju) = M - - [ du' p(u')\u-u'\ (22) 



J n x — 

P 2 2 J; 

The exchange term, X , explicitly mixes transverse and confined states, 

1/3 „ 

x[c a ^j a ,r s ] = -£^- / ^Ip^HI 2 ^^,^,^ 1 /^^) (23) 

where p ab {v) = J R du p a b(u) exp(-ivu) , p a b(u) = ^a{u)^pb{u), Go = V% (unpolarized) and G\ — 1 (polarized). The 
exchange function Y (see Appendix [A]) is given by 



Y(c a ,c b ,v) = ~ [ d 2 k f d 2 k' J 

77 J\k\*<c a J\k>\2<c b \k-k'\ l + v l 

and introduces a smooth variation in X as a function of r s . 

The Hartree-Fock ground state is determined by minimizing the total energy, Eq. (19), with respect to p, c a , and 
ip a , at fixed density, r s . We simplify this rather complex optimization problem, by considering only the completely 
polarized or unpolarized electron gas. For fixed concentrations, c a , the minimum of E with respect to ip a is independent 
from the in-plane kinetic energy. From the formal variation of the energy with respect to ip a we obtain 

— = — 3 H c a ^ a + -Y: V!f{u)i> h (24) 

^ a Ts S b 

H = -dZ + v p (u) (25) 
Vaf{u) = -~ / dvp ah (y)Y{c a ^r\l*G v vy™ (26) 

47T J R 

In the limit of small r s , the exchange energy is negligible, and the wave- functions, ip a = ip® are entirely determined 
by minimizing £, or, equivalently by the Hartree-equation 

Hoi>° a = (27) 



which leads to 



cO y0 

fa 



where X 7 ° n = X[c a , ip°, r s ] and 

S^=S[c a ,^° a } = -2j2c a K-v p (0)- fdup(u)v p (u) (29) 

is independent of r s . This provides us with a semi-analytical approximation for the total energy, E^, which appears 
to be very close to the full minimization of the energy including exchange, E m , for the densities considered. Whereas 
the in-plane kinetic energy term does not influence the shape of the distribution in z, it favors multi-mode occupation 
in the high density limit, r s — > 0. 



IV. SINGLE MODE SOLUTION 



For a single mode, we minimize Eq. ( 19 1 with p(u) 
equation, Eq. (24), 



(Ho 



i/jo(u) (cq = 1), so that tpo satisfies the non-linear Schrddinger 

(30) 



rl /3 V T)i>o = ~Wo 



and we obtain 



E x 



K 



Si 

4/3 



(31) 



where £\ and X\ are the values of the Hartree and exchange term using the optimal il>o- 

For the strictly two-dimensional electron gas, we have p^ 2D \u) — 5(u). Neglecting the zero-point energy of the 
confinement, E%, in this limit, and using f R dv Y(l, 1, v) — 32/3, we recover E^ 2D \p = 0) = l/(2r 2 ) — 8/(3-7rr s -\/2) for 
the unpolarized and E^ 2D \p = 1) = 1/r 2 — 8/(3nr s ) for the polarized electron gas. 
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A. Hartree solution without exchange, £*( 



We determine the one-mode solution of the Hartree-equation, Eq. (30) with Vqq c = 0, which determines the density 



distribution of the mode with Ao > in the high density region, r s —> 0. Assuming ipa(u) to be an even function of 
u, we restrict the discussion to u > in the following. From the leading order behavior A of the solution at large u 
where v p (u) vanishes, we use a series in Ao as ansatz for ip(u), 



^oW=a 2 ^(-l)^l„ 2H1 



(32) 
(33) 



fc>0 



where a = vAq, /q, and a k are to be determined (ao = 1). The density is then given by 



p(u) = a 4 *^2(-l) k p k Xo k+2 with p k = ^a.,a fc _j 



(34) 



k>a 



and the potential is obtained by integrating twice v''{u) = —p{u) for u > with the conditions v p (oo) = and 
<(oo) = 0: 



v p {u) = -a 2 ^2(-l) h v k X$ k+a with 



Vk 



Pk 



fe>0 



4(fc + l) 2 



and 



- a 4 £(-l) fc a fc (2fc + l) 2 A 2fc+1 

k 

k 

,(u)ip (u) = -a 4 ^(-l) fc w fc A 2fe+3 with w k = ^2v j a k - j 

fc>0 j=0 



(35) 

(36) 
(37) 



Thus, imposing %//q — (v p + a 2 )^j = leads to the equation: 



= a 4 A 



J2(-l) k a k m + If l]A 2fc + £(-l)*«, k X* 



2k+2 



k>0 



fc>0 



(38) 



From the definition ao = 1, we get po = 1, vq = \ and wq = \- The other terms are obtain by recurrence: 



a-k = 



4 • 

ttffc-1 



4k(k + l) 



(39) 



With these definitions all coefficients a k , pk, v k and w k are positive. The two parameters a and /o are determined by 
imposing ^o(O) = and the normalization: 



= ^(-l) fe (2fc + l)a fc / o fe cx^(O) 
) 



fe>0 



1 =n , 3 y (-l)Wo 

2 4^ 2(fc+l) 



fc>0 



(40) 
(41) 



The numerical results are given i n ta ble II together with the values of the different contributi ons to the Hartree- 



energy. In particular, from Eqs. (25 1, ([27 1, we have Aq 



((-A) + (v p )) and from Eq. (29), we have £f 



-2A - v p (0) - (v p ) with (v p ) = J M du p(u)v p (u) = -a 5 I] fe > (- 1 )' i; ' r fc/o + 2 ); T k = J2 J= o Pj v k-j, and the kinetic 

energy writes (-A) = (fa |-d 2 | Vo) = « 5 E fe >o r fe(" fo) k+1 1 '(* + 1) with r£ = £* =0 (2j + ^?^k- r Notice, that £° 
is independent of the polarization. 
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TABLE I. Parameters and various quantities of the single mode solution of the Hartree-equation without exchange term. 



fa 


15.5610024546998 


a 


0.465180466326271 


A 


0.216392866251527 


v(0) 


-0.674164469749883 


Mo) 


0.522553284700250 


(v) 


-0.307947186951202 


(-A) 


0.0915543206996701 


e? 


0.549325924198031 




FIG. 1. Comparison of the energies of the Q2DEG in the different phases within HF: black for the single mode unpolarized 
(p = 0) electron gas, red for the single mode polarized (p = 1) electron gas and blue for two occupied excited modes in z. 
For each phase we compare the HF energy using the Hartree density profile in z with the full HF minimization: dashed (resp. 

(resp. 



dotted) lines stand for E\ (resp. E\) of the single mode solution fro m Eq . 



(resp. 



E° 2 ) 



from Eq. 55 (resp 



43 



filled (resp. open) symbols stand for 
with squares (resp. diamonds) for the 



the energies including two occupied modes, Ei 
unpolarized (resp. polarized) gas. The red arrows indicate the transition between the unpolarized gas and the polarized gas 
at r a ~ 4.45 in the approximation using the Hartree density profile; minimization of the full HF energies shifts the transition 
to slightly higher density, r a ~ 4.05. Blue arrows indicate the transitions from the single mode system to two occupied excited 
modes increasing the density. The inset shows the transition region of the unpolarized gas. 



B. One-Mode exchange-energy in the Hartree-approximation, X° 



From the Fourier transform, p(v), of the ground state density, po = V'oi obtained from the Hartree equation, we 



23), to the total energy. Since we have Y > and < p(v) < 1, the 
gas is greater than its strictly two-dimensional value obtained with 



can estimate the exchange contribution, Eq 
exchange energy of the quasi-two-dimensiona 
pP D )(u) = l. 

The main contribution of the exchange-integral comes from the logarithmic singularity of the integrand at v = 0; 
details on the numerical evaluation are given in Appendix [B] and the results for the total energy are shown in FigjT] 
For densities corresponding to 0.5 < r s < 5, the exchange integral, X®, is well approximated by X®(p = l,r s ) w 
-0.4356 - 0.06127 ln(r s ) for the polarized gas. 

Within the Hartree-approximation, p(v) is independent of r s and polarization, p, so that a simple relation between 



o 



0.1 



1 t: ++++++++++++++ 



g -0.1 



o 



■0.2 



-+++++++ 



+ (£i-fi°)r s 2/3 




0.5 



FIG. 2. Importance of the the full minimization of all the different components of the HF energy, from Eqs. |29| and |43| for the 
polarized single mode gas. Shown are the total energy in Hartree times rf , together with the Hartree and exchange contributions 
in the same units. The gain in exchange energy is roughly twice the increase of the Hartree energy. 



X± of the polarized and unpolarized electron gas at different r s can be established 

*°(p=l,2>/§r,) 



#i(p = 0,r,) = 



V2 



(42) 



Using our approximate expression for X®(p — l,r s ) together with Eq. (42) in Eq. (31), we can estimate, that for 
r s > 4.56 the polarized phase is energetically favorable compared to the unpolarized phase. 



C. Full Minimization 



The full minimization assuming a single mode is done numerically (see Appendix [C] for the numerical details), and 
we have: 

+ ,43) 

where E\ and X\ depend on r s and the polarization, p. Figure [2] illustrates the small improvements due to the full 
minimization compared to the Hartree-approximation, E®. 



V. TWO-MODE SOLUTION 



In this section we look for the ground state energy with two modes, where the density is given by 

p(u) = (1 - c)Vg(u) + ci>l{u) (44) 

and c = ci is the concentration of the excited mode, ip\. Analogous to the discussion of the single-mode solution, we 
first minimize the Hartree-energy for given c with respect to ipo an d ipi to obtain £®- Then, we evaluated the exchange 
term within this solution, X® , and, finally, we minimize the full Hartree- Fock energy including the exchange. 



A. Two mode Hartree solution without exchange, 



Generalizing the single mode solution of the previous section, we express the wave functions as series of exponentials. 
Assuming ipo(v,) (resp. ipi(u)) to be an even (resp. odd) function of it, we restrict ipo and ipi to non-negative arguments 
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0.25 

0.20 

2 0.15 
u 

0.10 
0.05 



0.00 
0.00 



: . c(s) 




■ Fit: 1.537s 3 




- — lit: s°(l. 537 + 2.21s )/(l — 1.41s ) 






(a) : 



0.01 



0.02 



0.03 

o3 



0.04 



0.05 




FIG. 3. Graphical representation of the parameters for the two-mode model, see Eqs. 45|46 as functions of s = ^/Ai/Ao with 
Ai the eigenvalues of Ho, Eqs. 27 . (a): concentration c=c\ in the first excited mode versus s 3 . (b): fo and /i versus s. 



in the following 



A = Vfoe- au 



■00 (u) 

■01 (u) 



E v 2k+1 Y 2k' 



fc,fc'>0 



(45) 
(46) 

(47) 
(48) 



k,k'>0 



and 6 fc fc , are functions 



where a 2 = Ao, s 2 a 2 = Ai and a = 6 = 1. As shown in AppendixlDl the coefficients a k k 
of s only and can be determined by recurrence relations. Imposing theboundary conditions at u — 0: ip'o(0) — and 
01(0) = provide two equations independent of c and a: 



k t k' 



"k,k 



(2k + 1 + 2k' S )tifl 



k,k'>0 



(49) 



(50) 



fc,fc'>0 



In practice the series are restricted to k + k' < n. At large enough n, for fixed s, this system of the variables {fo, fi} 
has only one converging solution for positive fo and f\. The convergence with n depends on s. Relative convergence 
of one percent is reached at order n ~ 40. This slow convergence is due to the difficulty to fulfill the conditions at 
u = as we get close to the radius of convergence of these series. Machine precision is obtained using n ~ 120. 
Then the normalizations of ip and ipi lead to two simple equations determining a and c: 



E 



J0) fk+1 fk 
Vktk'JO Jl 



1 - c ^ 2(k + 1 + k's) 

k,k'>0 y ' 



s 4 a 3 \ - 



c ^ 2(k + (k' + l)s) 

k,k'>0 y x 1 ' 



(1) fkfk 
Pk.k'JQJl 



k fk' + l 



(51) 
(52) 



D2 



where p k a k , is defined in Eq 

The variations of a are essentially linear and given by a = 0.4608 + 0.44 c excepted at small c where we add the 
residual correction: 10 _3 (4.26 — 9.32c)/(l + 41.5c). The variations of the other parameters are given in Fig. [3] c is 
essentially proportional to s 3 , fo and f\ vary wi thin a factor of two. 



on the concentration c. (see Figj§left). 



Within the Hartree approximation, £i/ m Eq. ( 29 ) is still independent of r s and of the polarization, p, but depends 
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FIG. 4. Two-mode model results, (a): energy 82(c) versus c = Ci the concentration in the first excited state from Eq. (291 
compare d to the one-mode solution wit h £? = £3(0 = 0). (b): variations of the concentration c m i n , versus r s , that minimize 
E2 of Eq 54 (open symbols) or E2 of Eq 55 (full symbol). Squares (resp. diamond) stand for the unpolarized (resp. polarized) 



gas. 



B. Two-mode exchange term with the Hartree-approximation, X% 



The two-mode exchange term for two modes reads: 



X 2 (c,r s 



1/3 

ry 

4-7T 



dv 



a=0 



Paa{v)\ 2 C a Yi 



' r 1/3 r 



\~p m [v)\ 2 Y 2 (cyj*G p v 



(53) 



where cq = 1 — c and c\ = c. We refer to Appendix-JX] for the definition and evaluation of the exchange integrals Y\ 
and Y2 which have logarithmic singularities for small v, and to Appendix-|l3| for the evaluation of the exchange term. 

Using the Hartree approximation to determine the shape of the wave functions, the total two-mode energy is 
approximated by 



X p [(l-c) 2 + c 2 ] | f0( c ) *°(c,rj 



+ 



4/3 



(54) 



At fixed r s , a descent with respect to c allows us to determine the concentration c^ in which minimizes E^c). At 
small r s , a minimum c^ in 7^ is reached (see Figj4|right), and c° nin decreases as r s increases. The concentration in 



the excited mode vanishes at a critical value 



1.394(1) for the unpolarized gas. For the polarized g as, as 



increases, c = remains a local minimum. At r s — 2.775, the energy of the two-mode solution with cL, n ~ 0.015 
crosses the single mode energy. Thus, within this approximation, we find a first order transition for the polarized gas 
with a jump in the concentration (see FigSright). 



Full Minimization 



The complete minimization of the total energy with two modes containing kinetic, Hartree, and exchange energy, 
is done by first finding the ground state energy at fixed {r s , c} similar to the single mode case: 



Eo 



K p [(1 - c) 2 + c 2 ] £ 2 (c,r s ) X 2 (c,r s ) 



4/3 



(55) 



Then, at fixed r s , the minimum, c m i n (r s ), of the energy is found from a direct Newton-descent on c. The variations of 
Cmin(^s) are close to c^ in (r s ). They only differ significantly close to the transition. We find a transition at r SjC = 1.30(1) 
for the unpolarized gas and r s c = 2.50(2) for the polarized gas. In particular, no first order transition subsists for the 
polarized gas. The variations of the energy ^(cmin) versus r s are close to ^(Cmin) ( see Fig{lJ|. 
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D. Existence of three mode solutions 



We have further extended the method to study the occupation of three- modes. Unfortunately, the series used for 
the Hartree-approximation do not converge down to u — 0. Nevertheless, the solution can be found numerically, and 
we find the three-mode solution more stable for r s < 0.75 (resp. r s < 1.6) for the unpolarized (resp. polarized) gas. 
Since the exchange contribution becomes less important for smaller r s , we do not expect significant modifications 
from the full HF minimization. 

Approaching the high density limit, r s — ¥ 0, we expect an increasing number of occupied modes. For m modes, 
assuming c, = 1/m, the kinetic energy is K p /(m r 2 ), an d the dimensionless Hartree energy is a function of m only, 

£m({ c i = 1/™}) = F{ m )i as can be seen from Eqs (28 29|. Minimizing the total energy, E m ps K p /(mr 2 )+F(m)/r^ 3 , 
we can estimate the number of occupied modes in the high density limit: 

m 2 F'(m) = K p r- 2/3 (56) 

— 1/3 

Assuming a linear behavior of F for large m, the number of occupied modes diverges as r s as r s approaches zero. 



VI. CORRELATION ENERGY WITHIN THE LOCAL DENSITY APPROXIMATION 



Up to now, we have considered the total energy of the system within the Hartree-Fock approximation which 
neglects many-body correlation effects. Within density functional theory (DFT), the correlation energy per particle 
for m modes, C m , defined as the difference between the true total energy and the best Hartree-Fock solution, must be 
a functional of the electronic density onlyS. Using the local density approximation (LDAjP, we can write 



C m = ]^dup{u)4 u [ri u (u)] (57) 

where £c'°[ r s' D ] ^ s the correlation energy of the homogenous, three-dimensional electron gas at the (three-dimensional) 
density n 3D a 3 B — 3/(4ir[r 3D } 3 ) expressed in terms of the three-dimensional electron gas parameter r 3D . Using n 3D = 
a p(u)du/dz we get r 3D (u) = [3/8/d(u)] 1 ^ 3 r^ 9 . An estimation of the correlation effects is obtained by using the HF 
density, p(u), of the one and two mode density distribution, together with the Perdew-Zungei^ parametrization of 

<V y K n ■ 

Around the transition between one and two excited modes of the unpolarized gas, r s < 1.3, correlations, Eq. 
(57), lower the energy by typically less than 1%. Since the corresponding total density profiles (see Fig. [5| are 



smoothly varying with r s and with the concentration in the first excited state, c\, we do not expect important 
qualitative and quantitative modifications due to correlations in this density region. Energy minimizations including 
the LDA-correlation potential, V c [/>(it)] = SC m /8p(u), in the effective Schrodinger equation, confirm that Hartree-Fock 
accurately describes the high density region where the transition from single to two-mode occupation of excited modes 
occurs. 



VII. CONCLUSIONS 



We have studied the model of a quasi-two-dimensional electron gas where electrons are confined by a positive charged 
background localized in the plane z = 0. Similar to the 2DEG, the electronic density (r s ) is the only parameter of the 
system, however, the phase diagram is different due to possible transition from single to multi-mode occupation in 
z. Here, we have restricted the discussion to the most simple phases in the metallic regime neglecting the possibility 
of charge ordering and Wigner cry st allizat iorPi^J Already assuming a simple Fermi liquid wave function in the high 
density region, r s 0, we have shown that a transition from a single to two or more occupied modes in the confined 
direction takes place. Indeed, we expect that close to r s — three-dimensional features to be much more pronounced, 
as the dominant kinetic energy favors multi-mode occupations. Further, within HF, the transition between the 
polarized and unpolarized gas at r s ~ 4 occurs in between the correponding transitions of the 2DEG (r s ~ 2) and the 
3DEG (r s - 5 ) u . Similar to 2DEG and 3DEG, it is likely that the ferromagnetic phase of the Q2DEG is unstable 
against Wigner crys talliz ation within HF, however correlations are expected to stabilize the ferromagnetic fluid phase 
in higher dimension d 3 * 12 !, so that the spin-ordering of the Q2DEG may essentially differ from that of the 2DEG in the 
low density region. 

Within the Q2DEG, we expect that general aspects of the interplay between correlations and dimensionality can 
be studied without the need of a detailed microscopic modeling of a particular experimental device. This is of 
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FIG. 5. Comparison of the unpolarized charge density profiles p(u 

solutions, Eqs. (34 1 and (Dll while black lines stand for the optimized densities, as described in IV C and VC In 



rl /3 kFz). 



Red lines stand for the analytical Hartree 

(b) ci is 



the first-excited-mode concentration minimizing the total energy; at this scale, the corrections coming from LDA-correlation 
energy are negligible (b). Notice that the red line at r 3 — 1.3 in (b), with a rather small value of ci is close to the red line in 
(a). ' 



particular importance, since many experimental observations in quasi-two-dimensional electronic systems reflect strong 
correlation effect^. Up to no w prec ise calculations of correlation effects using quantum Monte Carlo methods have 
mostly be done for the 2DEGP 3 * 5 * 14 !, but perturbative inclusion of the underlying third dimension have shown to 
introduce important quantitative changes, e.g. concerning the spin susceptibility^. Within the Q2DEG model non- 
perturbative calculations are possible, and phases not contained in the 2DEG can be observed. As a side effect, a 
quantitative study of the Q2DEG using quantum Monte Carlo methods, may also provide a reference system, which 
is strongly inhomogeneous in one direction, so that, within DFT, corrections to the local density and generalized 
gradient approximations (GGA) should be more pronounced, and functionals beyond LDA/GGA can be tested (see 
refP). 



Appendix A: Properties of the exchange function Y 



The exchange function is given by the following integral: 

2 



Y(c a ,Cb,v) 



d z k \ d l k' 



1 



71-2 7|fc|2 <CQ J\k>\l<C b {k-k') 2 +V 2 

This function is positive for all v, even in is, and satisfies Y(c a ,Cb,v) = Y(cb,c a ,v) as well as aY (^r, 
Y (c a ,Cb,v). We find: 



(Al) 



Y(c a ,Cb,v) 



2lT 

dkk I dk'k' I d6 



1 



dkk / dk'k' 



k 2 + v 2 + k 12 - 2kk' cos(6>) 
1 



sj{k' 2 - k 2 + v 2 ) 2 + 4k 2 v 2 



dk k 



tanh 



c fa + v 2 - k 2 



_ x v 2 -k 2 



tanh 



^(k 2 - c b + v 2 ) 2 + <±c h v 2 k 2 + v 2 



= 2 dk 



,1,-1 c b + v 2 -k 1 v 2 

tanh — , in — 

^{k~c b + v 2 ) 2 + \c h v 2 2 k 

Y 2, , X -C a +C b + V 2 X + C a -C b + V 2 

X - c a — Cb — v + 2c a In —5 h 2c b In ■ 



2v 2 



2v 2 



(A2) 
(A3) 
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where 

X = +c a + c b ) 2 - 4c a c b (A4) 

In particular, within the context of the single mode solution it is convenient to introduce the function Y\{v) given by 

Yiiu/y/c) = Y(c, c, v)/c = Y (1, 1, is/Vc) 

-1 



with 



r 



11/4 



X + v 2 - 1 + 2c 



2c In 



X + v 2 + 1 - 2c 



YiO) = 2t-2-41n* 

whereas for the two mode model, we define Y% (y) 
Y 2 (c,v)=Y(l-c,c, v) 

Y 2 (c,v)=X-l-^ + 2(l-c)\n 2]/2 , _. 2/J 

Both, Fi and Y2, have a logarithmic singularity at v — and behave as v~ 2 at large i/: 
Yi{y) = -2 -41n|i/| +4|t/| + 0{v 2 ) 
%(c,v) = -4c In \v\ +2((1 -c)ln(l-c) - (1 - 2c) ln(l - 2c) - c) + 0{v 2 ) 

2c(l-c) c(l-c) 



+ 0(^- 6 ) 



(A5) 



(A6) 

(A7) 
(A8) 

(A9) 
(A10) 



Appendix B: Evaluation of the exchange term 



For the one- mode exchange term of the energy, X\ , we need to evaluate the following integral 



2?r J 



(Bl) 



where j3 



„V3 



and p(^) is the Fourier transform of p(u). For X®, the density is defined in Eq. (34) and p(v) can be 



computed from 



p{v) = 2 / dup(u)cos(vu) = 2a 4 YV-l) fe /o + Vfc 
Ja k>0 



2(fc + l)a 



4(fc + l) 2 a 2 + i^ 2 



To remove the logarithmic singularity of the integrand at v = 0, we introduce an auxiliary function e\ 



X x = - 



1 

2n 



ei (/3,r)+ / dv 




with 



ei(i/, P, t) = (e- Tl/ (l + rz/)) 2 /3(-41n(G p ^) - 2 + 4G p 0i/) 



ei (/3,r) 



dvei(v,P,T) 



9G P P 2 _,_/?/.. - , 
"IT 5 



57 



(B2) 

(B3) 

(B4) 
(B5) 



where 7 is the Euler constant and r = V 2 /° (2) is determined from p(z/) = l-pWv 2 + 0(v 4 ) {p {2) = 1.617362956587058 
using the solution of Eq. p4|). The integral in Eq. (B3) is then free of singularities and can be evaluated without 
major difficulties. 



Calculating the first integral in Eq. (53) contributing to the exchange term of two modes, X2, we adapt the above 



procedure for the integrals involving p aa using e%(v, r): 

e 2 (i/,/3,T,c) = (e- riy (l + ™)) 2 c/3(-41n(G p ^/v/c") -2 + 4G p Pv/y/B) 
e 2 (/3,r,c) = c 3/2 ei (/3,r^) 



(B6) 
(B7) 



The logarithmic singularity in the second contribution containing Y2 in Eq.(53l is cancelled by p a b{v) which is pro- 
portional to v at small v. 

Similar auxiliary functions are used to evaluate V^* c (w), whereas the logarithm singularity of Y%{v, c) in V°£ c (u) is 
again cancelled by p a b{v) oc v. 
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Appendix C: Details on the numerical minimization scheme 



Here we describe some details on the numerical minimization of the total Hartree-Fock energy, Eq.(19|. For 
simplicity, we restrict the discussion to the single mode solution where the fromal derivative is given by dip = Hd^ 
(Hd = 4"Ho/ r s^ 3 + 4Vbo/?' s )- We proceed using a quadratic minimization scheme. Let ip^ be the solution at step n 
and {dip^ -1 ^, dip^} the derivatives at step n — 1 and n. Energies E(e\,£2) are computed at ip^ +eidip( n ~ 1 ' 1 +£2^ 
for the six points (£1,^2) = (0,0), (±e, 0), (0, ±e) and (e, — e). By assuming a second order polynomial in e\ and £2, 
the minimum of E(£\ 1 E2) is determined analytically, and defines the solution at step n + 1. 

All functions of u, e.g. ip(u), are computed on a grid of 2 P points (i — io + 1)5 with i from to TP — 1, «o = 2 P ~ 1 and 
<5 = HmaxAo- Fast Fourier Transform (FFT) are used to compute p(v). In order to achieve good convergence small 
values of 5 are needed to accurately calculate the kinetic energy of the direct (Hartree) potential, whereas a small 
step in v is needed for the exchange energy which implies large values of u max - We found that u max = 150 and p = 10 
are good starting values at sufficiently large value of c. At small c, the spatial extension of the excited mode increases 
significantly which prevents accurate solutions for c < 10~ 3 . Interpolating ip(u) allows us to increase p at fixed u max - 

Appendix D: Recurrence relation for the two mode Hartree solution 

We determine the recurrence relation of the series coefficients in the two mode case. The densities are given by 

E Y 2k Y 2k' n _ JO) , 4 .(1) / n1 \ 

Pk,k>X X 1 p k . k , - p k _ lk , + S p kk ,_ x (Dl) 

k,k'>0 

Pk.k' 



4 

p = a 



^2 Yl I' Phi' = 53 b M' b k-j,k'-j' (D2) 



J=0j'=0 j=0j'=0 

with the convention that a^i k > — b k _i — 0, and the potential is defined as 



v p (u) = -a 2 VkvXl k xf with v ktk , = 7 fc , (D3) 



We have 



k k' 



v p ip = — y= w kM ,Xl k+1 Xl k ' with w k ,k' = v 3,i ,a k-]M'-r (D4) 

V C k,k'>0 j=0 j'=Q 



v p ^ = ~ £ <,k>Xl k Xl k ' +1 with "V.r ; > V (D5) 



k k' 



IC 

k,k'>0 j=0 j'=0 



Imposing ipQ — (v p + a 2 )V>o = and ip'{ — (v p + sa 2 )ipi = gives: 

]T a k>k , [(2k + 1 + 2fc'.s) 2 - 1] X^xf + w k ^X^xf (D6) 

k : k'>0 k,k'>0 

b k ,k> [(2fc + (2k' + 1) S ) 2 - s 2 ] A 2fc Af'+ 1 + w'kM'X 2 o k xf +1 (D7) 

k,k'>0 k,k'>0 



with the following solution for (k,k') 7^ (0,0): 



w k,k' (T)fl , 

ak > k ' ~ A(k + k's)(k + l + k's) [ ' 



IV 



h ' k '~ 4(k + k's)(k'+(k' + l)s) (D9) 
Thus, the coefficients a k . k > and b k> w, as well as p k . k > and v kyk i, are rational functions of s only. 
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